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1 Introduction 


The theory governing the strong force, QCD, predicts distinctive changes in the behaviour 
of hadronic matter when exposed to extreme conditions, i.e. at high temperature and/or 
at high baryon density. It is the aim of large experimental programs to create matter in 
such an environment and to study its properties. In order to interpret the experimental 
hndings, one of the important questions to be answered is to quantify how properties of 
hadrons produced in heavy ion collisions are changed in a hot medium. 

In this paper we investigate possible changes in the internal structure of hadrons at a 
temperature below the critical temperature T^, i.e. we look for precursors of the phase 
transition for the pion. More specihcally, we use lattice QCD to calculate spatial two- and 
three-point functions for a pion, to extract the pion ’form factor’ at hnite temperature. 

Finite temperature form factors of the pion have been considered in the framework 
of a variety of theoretical approaches. An increase in the charge radius with increasing 
temperature was found in the Nambu-Jona Lasinio model pp. Dominguez et al. j21 used a 
hnite energy QCD sum rule to show that the charge radius of the pion increases with T and 
diverges at some critical temperature. An extensive study of the pion electromagnetic form 
factors at hnite T was undertaken by Song and Koch j3] in the time-like region. Using an 
ehective chiral Lagrangian to one loop, the form factor at T > 0 was found to be reduced 
in magnitude, in particular in the vector meson dominance region. Similar results were 
recently obtained by Nicola et al. |3] who work in chiral perturbation theory to one loop. 
They calculated the form factor in the space-like region and found that the charge radius 
initially stays constant, but then increases at higher temperature. 

There have already been several lattice QCD investigations of some properties of a pion 
embedded in a heat bath. A recent review of these studies, both quenched and unquenched, 
can be found in Ref. j3]. Most of the work considered spatial correlators and extracted 
the so called ’screening masses’ and spatial ’Bethe-Salpeter wave functions’ of the pion. 
The general conclusion was that for temperatures below the phase transition temperature, 
Tc, there are no signihcant differences between the finite temperature ’screening’ mass 
and the free pion ’pole’ mass, extracted at T = 0. Furthermore, the ’Bethe-Salpeter wave 
functions’ extracted from spatial two-point correlators at T < Tc were very similar to those 
at T = 0. However, as subsequent calculations Pd of the pion form factor at T = 0 have 
shown, conclusions about the internal structure based on these wave functions obtained 
from two-point correlators are not very reliable. 

We therefore study an additional observable, the three-point function of the electromag¬ 
netic current between pion states. This enables us to extract the electromagnetic vertex of 
the pion at T > 0 and to draw direct conclusions about the internal structure. We do this 
for a medium which has vanishing net baryon density and a temperature of T = 0.93 Tc. 
Our study of the pion vertex function for space-like photons is the hrst such investigation 
of the spatial structure of a hadron at hnite temperature with lattice QCD. For our study 
we use an 0(a) improved Wilson action, together with the consistently improved vector 
current in order to ensure the absence of 0{a) effects in the matrix element. 

In Sec. II, we hrst discuss some general features of the meson form factors at hnite 


2 


temperature followed by some technical details in Sec. III. In Sec. IV, we present the 
calculation of the pion two-point function and discuss our results for screening mass and 
dispersion relation. The pion electromagnetic vertex is then extracted from the three-point 
function and discussed in detail in Sec. V. Sec. VI contains our summary. 

2 Formalism at T > 0 

Temperature is introduced in the path integral formalism by restricting the Euclidean time 
direction and imposing (anti-) periodic boundary conditions. In the lattice approach, the 
temperature is then dehned through T = {N^- a)~^. To facilitate the determination of the 
pion form factor on the lattice, one has to calculate two observables. These are the two- 
and three-point Green’s function of an interacting quark-antiquark pair, which carries the 
quantum numbers of a pion, at large separation. Since at higher temperatures the temporal 
direction is rather short, it becomes difficult to reliably hlter out the ground state from 
correlators in the r-direction and we use instead spatial correlators in the z-direction. 

2.1 The two-point function 

The two point function is given by 

G{x, t) = {(j){x, r)0^(O, 0)) , (1) 

where for dehniteness we assume that we are dealing with a tt’*' meson, 

(j)\x) = ^pu{x) 7 ® -ipdix) . (2) 

In the following, we will suppress all flavor, SU{3) color and spin indices. We want to 
study the spatial correlator in the 2 ;-direction, 

G{z,p) = J dr J dx J dy G{x,t) , (3) 

where 

X={x,y,T) P={Px,Py,P4) ( 4 ) 

denote the three-vector parts of the coordinates and momenta in the so called ’funny space’ 
[Hj. Due to the periodic boundary conditions for the pion, p 4 is restricted to the Matsubara 
frequencies, 

P4 = 27inT = UJn ■ (5) 

In the following we consider only the lowest contribution, (Uq = 0, since the next mode at 
uji ~ 1.6 GeV is already quite heavy. 
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2.1.1 Dispersion relation and wave function renormalisation 

The inverse pion propagator in Euclidean momentum space can be written as 

A-\p,P4;T) , ( 6 ) 

where m is the bare pion mass and effects due to the presence of the heat bath are incor¬ 
porated into the self energy 11. In contrast to the situation at T = 0, the propagator can 
now depend separately on 


p'^ = p]_+Pz and Pi = n ■ p. (7) 

Here, the four-velocity of the heat bath, is given by 

n^ = (0,0,0,l). (8) 

In terms of the momentum space propagator in Eq. the spatial correlator in the 2 ;- 
direction is given as 

G(z,p)^ I A{p,pt-,T). (9) 

Its behavior at large 2 ; is determined by the poles with the lowest pi value of the propagator, 
Eq. (|n)). For pi = 0, the poles occur when the spatial momentum satishes 

p2 = _^2 _ 0; T) = -ml^ T , (10) 

where msc,T denotes the temperature dependent screening mass. For a given value of the 
transverse momentum pj_, the pole in p^ is therefore located at pz^, with 

pIo = - EUpL0;T) , ( 11 ) 

where we have introduced the screening energy, 

Escipl; T) = ^JplTm^ . (12) 

The state we hlter out for large separation 2 ; at a given transverse momentum pj^ and 
P 4 = 0 is thus the state with the lowest screening energy, which satishes the dispersion 
relation Eq. m or Eq. m- This state is in principle different from the ground state of 
the pion with the lowest energy. For simplicity, however, we will refer to the state with the 
lowest screening energy as the ’ground state’. 

In order to obtain the wave function renormalisation for this state, we expand the 
propagator around the pole p^ = Using 

n(p^0;T) = n(-m^^^r,0;T) 

+ (p^ + "^L,t) ^n(p^0;T)|p2=_^2^^^-F... , (13) 
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( 14 ) 


and Eq. m, one can write the inverse propagator as 
A-i = (p2 + + -^n{p\ 0; T) 

where the remainder vanishes at the pole, 


= 0 


Near the pole, we therefore obtain 


A(p,P4 = 0;T) = „ —^ , 

Pi + Pi + 

where the renormalisation constant Z is dehned by 

dp'^ 


Z ^ — 1 + 7——n(p^, 0; T)|p2=_m2 


sc,T 


(15) 


(16) 


(17) 


2.2 The three-point function 

The second observable we nse is the three-point fnnction. For the pion, nsing degenerate 
qnark masses, only connected diagrams contribnte. The qnark-antiqnark pair propagates 
from the sonrce at Zi = 0 to the sink at Zf. The photon conples to the propagating 
quarks at an intermediate point z. In the continuum, this function reads in the notation 
introduced in Eq. 

r^T i^i/T 

G^{zf,z,Pf,Pi) = / d^Xf / 

X {(l>Rixf,Zf)jf,{x,z)(j)^{0,0)), (18) 

where is the quark vector current to which the photon couples. Details of the current 
will be given in Section 3. 


2.2.1 Current structure and form factors 

Given the four-vectors n^pi and p/, a matrix element of the electromagnetic current oper¬ 
ator has for the pion the general Lorentz structure, 

Jf. = (7r(P/)|jpk(Pi)) 

= e^{{pi+pf)f, F + G + H} , (19) 

where 

= (P/ - Pi)^^ (20) 

is the photon four momentum and Ctt the pion charge. The functions F, G, and H are form 
factors which are functions of the independent scalar variables and are subject to conditions 
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arising from current conservation or the Ward identity. In contrast to the situation at T = 0 
they can depend on more scalar variables, namely 

and n-pf . ( 21 ) 

The last two scalars amount to and p/^ 4 . When calculating the three-point function, 
Eq. (fTRjl . we choose to project out states with transverse momenta of same magnitude, 

pi,/ = P\,i = pI • (22) 

By choosing large spatial separations of the vertex from both source and sink and P 4 = 0, 
we hlter out the initial and hnal pion in the ’ground state’ which satishes the dispersion 
relation Eq. Ha. It is therefore easily seen that the form factors with this choice of 
kinematics will only depend on two scalars, the screening mass msc,T and 

Q" = -<f = {P^j-P^,f- (23) 

In our applications the momentum transfer to the pion, is varied by changing the angle 
between p_|_ ^ and p_|_j and by choosing different values for Eq. (j22I)- 

Current conservation for matrix elements of the current operator between initial and 
hnal pion states yields in general 

0 = (P/ - pI) Q^) + + n ■ {pf - pi) HQ^) . (24) 

Since for our symmetrical kinematics p‘j = pf and n-{pf—pi) = 0, we see that G must vanish. 
Furthermore, the term involving H will not contribute if we consider spatial components 
of the current. Choosing in particular the 2 ;-component, we obtain for the lattice version 
of the pion current 

{'^{p±j) \jz\ '^{p±,i))lattice = e^F{Q^,ml ^ . 2 N 

2 VPzjiPl) PzApI) 

= A ■ (25) 

As already pointed out in Ref. [Zj, the cancellation of the kinematical factor in the last 
step due to our symmetric choice of momenta makes the extraction of the form factor F 
from the lattice data more reliable. 

2.2.2 Effective charge 

In comparison to the free pion current at T = 0, there are two modihcations at T > 0. 
First, the overall renormalisation constant Z, Eq. (ED at hnite temperature, might differ 
from the corresponding value at T = 0. Secondly, the vertex operator gets modihed. At 
the photon point, where = 0 and pj_j = Px j, we absorb both effects into an effective 
charge e^ff, 

AiPL)\jz\AP±)) = ee// 2 p, . (26) 
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The modification of the vertex can in this case be obtained from the Ward-Takahashi 
identity for the vertex operator 

{p' -P)ti = e^{A"^(p',P4;T) - A“^(p,P4;T)} (27) 

Taking hrst p' = p and then the limit —>■ Pz, one obtains at the photon point 

jz{Q^ = 0) = ^A"^(p,P4,T) 

OPz 

= 2e^pz (^l + ^n(p2,0;T)^ (28) 

When evaluating the matrix element of the vertex operator between states satisfying 
the dispersion relation Eq. we have 


jziQ" = 0) 


2 Pz 



d 

dp'^ 


n(p2,0;T) 



2 

sc,T 


= 2e^PzZ-\ (29) 

where Z was dehned in Eq. m- The vertex correction will thus precisely cancel the wave 
function renormalisation factor Z, resulting in 


^e// ~ • (20) 

Thus if we use pion states obtained from spatial correlators, the effective charge determined 
from spatial components of the conserved current is the same as the free charge. This 
complements the remark in Ref. [H], where it was noticed that the effective charge obtained 
from the 2 ;-component is different from the one obtained from the t-component if the same 
wave function renormalisation constant is used in both cases. This difference is due to the 
fact that Lorentz invariance is broken and the self energy depends separately on p and p 4 
at T > 0. 


3 Technical Details 

Most practical aspects of our lattice calculations are identical to the procedures in Ref. [7] . 
We therefore only give a short summary. 

Our calculations were done on a Nfj x Nr = 32^ x 8 lattice. The Nr value corresponds 
to T = 0.93 Tc at the chosen value of /? = 6.0. In comparison to Ref. [7j, the spatial size 
of the lattice has been increased from N„ = 24 at T = 0 to = 32 to have the same 
length in the correlation direction. We generated 0(200) conhgurations, twice as many as 
at T = 0, since we expect an increase in fluctuations in the vicinity of the phase transition. 
To facilitate comparison with the results in [7], we used the same hve k, values, 

K = 0.13230, 0.13330, 0.13380, 0.13430, 0.13480 (31) 
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which correspond to pion masses of 360 — 970 MeV^ at T = 0. The action we use is 
the Sheikholeslami-Wohlert action ng with the non-perturbatively determined m value 
csw = 1.76923, which is exact to order 0{a). 

In order to obtain matrix elements of the current that are correct to order 0{a), we 


have to use the appropriate vector current for the chosen action, 
current is [1210 

This improved vector 

ji = EvVi + acy , 

(32) 

where 


Jm = 7m 

(33) 

is the local quark current and 


Ti,u = i’ix) ^^{x) . 

(34) 


The renormalisation coefficient Zy and the improvement constant cy have been determined 
non-perturbatively ug. This improved current guarantees that on-shell matrix elements 
only have 0{a'^) errors. For comparison, we also use the conserved lattice current in 
our calculations 


(1 - 7m) U^{x) ^{x + /i) - ^{x + /i) (1 + 7^) ul{x) ^{x) , (35) 

which has 0(a) discretisation errors away from the forward direction. To enhance the 
contribution from the pion ground state at the sink point {xf,Zf), we use an extended 
operator 0^ with a suitably chosen inter-quark distance R. To keep the calculation gauge 
invariant, the separated quarks at the sink are connected by gauge links which are ’fuzzed’ 
to simulate the tube-like nature of the gluon cloud. In these steps we follow the scheme 
developed by Gupta et al. UTI and Lacock et al. the fuzzed gluon links at the pion sink 
are created with a link/staple mixing of 2 and a fuzzing level of 4. As a measure for the 
effectiveness of this method, we used the speed with which the effective screening energy. 


EfMz,P) = In 


{Gr{z,p)) 
{Gr{z -I- l,p)) 


(36) 


stabilizes for increasing We varied the quark separation R and found that, just as at 
T = 0, the value R = 3 was optimal. The same R was also used for the three-point 
functions. 


4 Results for the two-point function 


In order to extract information from the lattice data, the two-point function, Eq. ©, is 
parameterised, with p 4 = 0 as 


1 

Gr{z,p^) = \/zS(pi) Zo"(Pi) e-®-" ‘ol> cosh 


n=0 


N 


^We use a = 0.105 fm from to set the scale. 


( 37 ) 
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Table 1: Fit results at T = 0.93 compared to corresponding data at T = 0. 


K, 

msc,T 

{T = 0) 

K 

mv 

mv {T = 0) 

mp uni 

0.13230 

0.511(3) 

0.516(2) 

0.13230 

0.574(11) 

0.587(19) 

0.623(2) 

0.13330 

0.410(4) 

0.414(2) 

0.13330 

0.520(13) 

0.528(17) 

0.550(2) 

0.13380 

0.353(4) 

0.356(2) 

0.13380 

0.495(15) 

0.501(19) 

0.515(3) 

0.13430 

0.283(5) 

0.287(3) 

0.13430 

0.466(18) 

0.477(21) 

0.485(3) 

0.13480 

0.187(8) 

0.194(4) 

0.13480 

0.431(25) 

0.454(43) 

0.448(13) 


(a) Pion masses. (b) VMD fit parameter mv compared to the free p- 

mass. 


where = 32 is the extension of the lattice in the 2 ;-direction. The cosh-structure 
is due to the periodic boundary conditions in 2 ;. The quantities denote the matrix 
elements n, p_|_)p , the overlap between the trial state and the pion state 

with screening energy E^^{pY) and momentum p^. Since we use smearing techniques at the 
sink to suppress higher excited states, the parameterisation can be restricted to two states. 
For a given momentum pj^, the parameters we £t for each state are thus the screening 
energy, and the product of the amplitudes, Zq, in Eq. (jSID. For the momenta, we 
choose 

pI = 0 , 0.039, 0.077, 0.154, 0.193 , (38) 

which are integer multiples of where Pmin = These values include the momenta 
used for our form factor calculation. For each nonzero value of we use several pj^ 
with different directions in the {x, y} plane. On the lattice there is still some rotational 
symmetry left and the two-point function is independent of the direction of the transverse 
momentum. We therefore average over all momenta with the same length to increase the 
stability of the data. 

There might still be correlations between different conhgurations. To deal with these, 
we obtained all our results using the jackknife method 1201 in]. 

4.1 Screening masses and dispersion relation 

The results for the screening masses, msc,T-, obtained for pj^ = 0, are given in Table mg 
Also shown in the table are the free pion masses, mjr, obtained for the same action at 
T = 0 in Ref. [7j. As can be seen, the screening masses agree with the corresponding pole 
masses at zero temperature within error bars. This conhrms earlier work 12211201, where no 
signihcant difference between screening and zero temperature pole masses were found for 
a TT meson at temperatures below Tp, similar observations were also made for the p meson. 

In order to investigate the dispersion relation for the lattice states, we plot the screening 
energies as function of the transverse momentum in Fig. ^ The screening energies obey 
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Figure 1: i?sc(pl) for different pion masses; lines: continuum dispersion relation Eq. () 12 j) . 


the continuum dispersion relation, Eq. m, quite well. This nicely conhrms that lattice 
artifacts are largely suppressed since we are working with an improved action, which only 
receives corrections in The error in Esc{p±) grows with increasing momentum and 

decreasing pion mass. Due to large fluctuations, reliable results for the lowest pion mass 
could not be obtained at the two highest momenta. 

5 The three-point function and the vertex function F 

For the three-point function, we consider a pseudo-scalar source at 2:4 = 0, a sink at 2 ;/, 
and a coupling of the photon at 0 < 2 ; < 2 ;/. As was discussed in Ref. [Zj, for T = 0, the 
most reliable method to extract the form factor is the use of simultaneous hts of the two- 
and three-point functions. Consequently, this will also be the method used in this case. In 
analogy to the zero temperature analysis, we have varied the £t range in both correlation 
functions to investigate systematic uncertainties. We found no signihcant changes in the 
ground state parameters. 

We carried out our simulations for the hve k values in Eq. m, corresponding thus 
to the different screening masses in Table |2(a)| Furthermore, we chose three external 
momenta, p]_ = 0.039, 0.077 and 0.0193. As can be seen from our discussion in Sec. 12 . 2 . II 
we do not expect a dependence on for our specihcally chosen kinematical situation. 
The data sets for the different transverse momenta will therefore yield an indication of 
numerical instabilities in our results. For the lowest pion mass and highest p\, fluctuations 
overwhelmed the data and we did not extract the parameters of the three-point function. 
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Figure 2: Renormalisation constants Zy as function of quark mass for different T and p±. 

As shown above, for our kinematics and by considering the 2 ;-component of the current, 
the resulting matrix element takes on a very simple structure with only one single vertex 
function F, which depends on and the screening mass. In our calculations, we have 
considered three choices of the current operator on the lattice: the conserved current, 
Eq. (j35|l . the renormalised local current, Eq. O, and the improved current, Eq. (E2D- 
At = 0, for periodic boundary conditions, the three-point function calculated with a 
conserved current should satisfy the relation pg 

G3(zf,z;p,p) -G3{zf,z',p,p) _ , 

G(zf,p) 

where Zf < z' < N„. In this form, the current operator is inserted twice to measure two 
separate contributions to the total charge reaching the sink at Zf. The first takes into 
account the charge from the source reaching the sink by passing through z. The ’second 
insertion’ at z' accounts for the charge that leaves the source in the negative z-direction and 
arrives at the sink by passing through z', which is possible due to the periodic boundary 
conditions. We have confirmed that for the conserved current Eq. (P|) is satisfied to high 
accuracy, typically to 0(10“^). 

For both the local and the improved current we need the overall renormalisation con¬ 
stant Zy, 

Zy = Z°(l -h byrng). (40) 

We have determined it by demanding that Eq. dsni) is satisfied. Since the additional tensor 
term in the improved current is a total divergence, it does not contribute to the total 
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Figure 3: The inverse of the vertex function F as a function of Q^. 


charge and therefore Zy is identical for the local and improved current. Fig. |21 shows 
the Zy-values we obtain for several k and p\ values. The higher momenta are more 
problematic for the extraction of the vertex functions. They have larger error bars and the 
extracted Zy values tend to he slightly higher than the lower momentum values. As the 
figure shows, the T = 0.93 results deviate by less than 2% from T = 0 values extracted 
in Ref. [7] . They also have the same linear dependence on the quark mass. Thus, our values 
for Zy as well as by are in very good agreement with the non-perturbative determination 
of Bhattacharya et al. us. In the following, we therefore use their value for the small 
improvement constant cy for the vector current, which contributes to the vertex only away 
from the forward direction. 

5.1 Results for the vertex function F 

We now discuss the values for the vertex function we extracted with the improved current 
for different photon momenta Q^. At T = 0, it was found that a monopole parameterisation 
of the pion charge form factor described the obtained results quite well for the lower 
values, a feature which can be explained by the vector meson dominance (VMD) model. 
We investigated if the analogous parameterisation, 

ml t) = [1 + ^]-\ (41) 

where my is a fit parameter, also yields a useful description for the vertex function we 
extract at T > 0. We therefore plot 1/F, which should be a straight line for the parame¬ 
terisation to work. 

As can be seen in Fig. 13 for two values of the pion (screening) mass, our numerical 
results show that for up to 0.4 a monopole £t will work quite well. At higher Q^, the 
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Figure 4: The vertex function for different pion masses. Curves: fits to the VMD model. 


lighter mass shows a stronger deviation from a straight line, but with larger error bars as 
well. When comparing data from the different transverse momenta sets, we see that the 
higher p]_ values have larger error bars and for the lighter mass are more scattered around 
a straight line. 

Fig.lUshows all our results, appropriately averaged over the different values, together 
with the monopole fits. For the lightest pion mass, the fluctuations for the higher p\ values 
were too strong to allow a reliable extraction of the vertex function at large values. As 
expected, F^Q"^) drops off more rapidly as the pion mass decreases. The fits deviate more 
strongly from the data at higher for the lighter pions. As shown in Table |2(b)[ the 
extracted parameter mv agrees quite well with the corresponding fit parameter extracted 
at T = 0. Both lie close to the free p-mass obtained from lattice QCD in Ref. This 
agreement with rrip gets better as the pion mass decreases towards the physical value. Our 
£t at T = 0.93 Tc thus in general supports the simple VMD picture for the low data also 
at a temperature just below T^. A more detailed look at the function F is taken in Fig. El 
The two different pion masses show how the difficulties to extract information about F 
increase as the pion mass decreases. However, within the error bars, both examples again 
show that there is no overall significant difference to the form factor F{Q‘^) we extracted at 
T = 0. This is in contrast to most theoretical expectations based on effective models. In 
terms of the mass parameter mv of our monopole fit, our results therefore do not support 
a signihcant dropping of the vector meson mass as T increases to 0.93 at least not for 
the state with the lowest screening mass we project out by our spatial correlators. 

By identifying the vertex function F{Q‘^, msc,T) aXT = 0.93 Tc with the form factor of a 
pion embedded in a heat bath, we can also translate our findings into statements about its 
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Figure 6: The pion ’radius’ at 0.93 T^. as a function of its mass. 


spatial extension and mean square radius. Fig. IHl shows (r^) as obtained from the slope of 
F a.t = 0, using the VMD parameterisation. Both the T = 0 and T = 0.93 Tc data lead 
to the same picture. There is no signihcant difference between the value for (r^) and in the 
dependence on the pion mass for both temperatures. Using the VMD parametrisation and 
extrapolating the ’p’ mass my, linearly in ml, to the physical limit, we arrive at a value 
for (r^) that is about 15% below the experimental value for a free pion [2^]. This small 
discrepancy could be due to the quenched approximation. However, Alexandrou et al. I2SI 
have calculated density-density correlations in quenched as well as unquenched QCD at 
T = 0 and found only rather small effects for the pion. Recent unquenched calculations of 
the free pion form factor itself can be found in mM 


6 Summary 

In this paper, we have presented the hrst investigation of the electromagnetic vertex of a 
pion at hnite temperature by means of numerical simulation. The temperature was chosen 
as 0.93 Tc, just below the critical temperature for the transition to the deconhned phase 
of QCD. 

We considered two- and three-point functions for the pion for various transverse mo¬ 
menta. As the temporal extent of the lattice is limited by the inverse temperature we have 
calculated spatial correlators in the 2 ;-direction to study the dispersion relation and to ex¬ 
tract the vertex functions of a pion embedded in a heat bath. For several quark masses we 
conhrmed earlier Endings that within errors the screening mass at 0.93 Tc is the same as 
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the free pion mass calculated on the lattice. Furthermore, except for the lightest pion and 
the largest momentum where errors are large, we could show that the dispersion relation 
for the lattice results was sufficiently well described by the continuum expression. 

For the symmetrical kinematical conditions specihcally chosen by us, the electromag¬ 
netic vertex of the pion can be described in terms of only a single vertex function nisc), 

which depends on the square of the photon four-momentum and the screening mass at the 
chosen temperature. For T = 0, this function is equal to the free pion form factor F(Q^). 
Within error bars, we found that for all quark masses and all we considered, the vertex 
function agreed with the corresponding form factor at T = 0. Furthermore, we showed 
that at low the vertex function for T = 0.93 Tc is well described by a monopole £t 
as suggested by a vector meson dominance model. The mass parameter my we extracted 
equals the fit parameter found earlier at T = 0 within errors. As the quark mass was de¬ 
creased, this parameter was seen to approach the lattice value rup for a p-meson at T = 0. 
The simple vector meson dominance model thus remains a good description of the low 
results also at a temperature close to T^. 

By interpreting the vertex function as the form factor of the pion state with the lowest 
screening mass, we can also translate our low results into a mean square radius < >. 

Contrary to most predictions based on effective models that the radius of a pion will 
gradually grow as the temperature of the medium increases, we did not find significant 
differences from the situation at T = 0. 

To conclude, our results for the screening masses at T = 0.93 Tc confirm what earlier 
studies of the pion had already indicated. Our main new finding is that there is also no 
significant change when one looks at the vertex function F, a quantity that is directly 
sensitive to the internal structure of the pion. It will of course be interesting to continue 
these studies of the three-point functions to higher temperatures to see whether significant 
changes occur when Tc is crossed. Such studies are presently in progress 


References 

[1] H. J. Schulze, J. Phys. G20, 531 (1994). 

[2] C. A. Dominguez, M. S. Fetea, and M. Loewe, Phys. Lett. B387, 151 (1996), hep- 
ph/9608396. 

[3] C. Song and V. Koch, Phys. Rev. C54, 3218 (1996), nucl-th/9608010. 

[4] A. Gomez Nicola, F. J. LLanes-Estrada, and J. R. Pelaez, Phys. Lett. B606, 351 
(2005), hep-ph/0405273. 

[5] F. Karsch and E. Laermann, (2003), hep-lat/0305025. Published in *Quark gluon 
plasma III, Hwa, R.C. (ed.) et al* 1-59. 

[6] J. van der Heide, M. Lutterot, J. H. Koch, and E. Laermann, Phys. Lett. B566, 131 
(2003), hep-lat/0303006. 


16 


[71 J. van der Heide, J. H. Koch, and E. Laermann, Phys. Rev. D69, 094511 (2004), 
hep-lat/0312023. 

[8] V. Koch, E. V. Shuryak, G. E. Brown, and A. D. Jackson, Phys. Rev. D46, 3169 
(1992), hep-ph/9204236. 

[9] R. G. Edwards, U. M. Heller, and T. R. Klassen, Nucl. Phys. B517, 377 (1998), 
hep-lat/9711003. 

[10] B. Sheikholeslami and R. Wohlert, Nucl. Phys. B259, 572 (1985). 

[11] M. Liischer, S. Sint, R. Sommer, P. Weisz, and U. Wolff, Nucl. Phys. B491, 323 
(1997), hep-lat/9609035. 

[12] G. Martinelli, G. T. Sachrajda, and A. Vladikas, Nucl. Phys. B358, 212 (1991). 

[13] M. Liischer, S. Sint, R. Sommer, and H. Wittig, Nucl. Phys. B491, 344 (1997), 
hep-lat/9611015. 

[14] M. Guagnelli and R. Sommer, Nucl. Phys. Proc. Suppl. 63, 886 (1998), hep- 
lat/9709088. 

[15] T. Bhattacharya, R. Gupta, W.-J. Lee, and S. R. Sharpe, Phys. Rev. D63, 074505 
(2001), hep-lat/0009038. 

[16] L. H. Karsten and J. Smit, Nucl. Phys. B183, 103 (1981). 

[17] R. Gupta, D. Daniel, and J. Grandy, Phys. Rev. D48, 3330 (1993), hep-lat/9304009. 

[18] UKQGD, P. Lacock, A. McKerrell, G. Michael, 1. M. Stopher, and P. W. Stephenson, 
Phys. Rev. D51, 6403 (1995), hep-lat/9412079. 

[19] M. Gockeler et al, Phys. Rev. D57, 5562 (1998), hep-lat/9707021. 

[20] M. H. Quenouille, J. Roy. Statist. Soc. B Bll, 18 (1949). 

[21] J. W. Tukey, Ann. Math. Stat. 29, 614 (1958). 

[22] E. Laermann and P. Schmidt, Eur. Phys. J. C20, 541 (2001), hep-lat/0103037. 

[23] E. Laermann et al, Proceedings of SEWM04 , 206. 

[24] K. Barad, M. Ogilvie, and G. Rebbi, Phys. Lett. B143, 222 (1984). 

[25] NA7, S. R. Amendolia et al, Nucl. Phys. B277, 168 (1986). 

[26] G. Alexandrou, P. de Forcrand, and A. Tsapalis, Phys. Rev. D66, 094503 (2002), 
hep-lat/0206026. 


17 



[27] D. Brommel et ai, PoS LAT2005, 360 (2005), hep-lat/0509133. 

[28] JLQCD, S. Hashimoto et al, PoS LAT2005, 336 (2005), hep-lat/0510085. 

[29] J. van der Heide, J. H. Koch, and E. Laermann, In preparation. 


18 



